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Abstract 



We propose a new approach to the mathematical description of hght propagation in 
a single-mode fiber light-guide (SMFLG) with random inhomogeneities. We investi- 
' gate statistics of complex amplitudes of the electric field of light wave by methods of 

, the random group theory. We have analyzed the behavior of the coherence matrix 

On ' of a monochromatic light wave and the polarization degree of a nonmonochromatic 

radiation in SMFLG with random inhomogeneities as the length of the fiber tends to 
infinity; in particular, we prove that limit polarization degree is equal to zero. 



fj^: 1 Introduction 

"j^ I To be able to analyze statistical parameters of light propagation in a single-mode fiber 

light-guide (SMFLG) with random inhomogeneities is vital in connection with broad appli- 
cation of SMFLG in promising directions of industrial development such as optical commu- 
nication iQ, ^ and improvement of sensors of various physical parameters |^ ^ ^, 0, ^ . 
^ I Existence of random inhomogeneities in real optical fibers leads to appearance of cou- 

■ pling of orthogonal polarization modes, it is accompanied by exchange of energy between 

modes. In turn, the mode coupling causes random variations of a polarization degree 
of radiation in SMFLG and, as a result, undesirable phenomena in devices and commu- 
nication lines developed on the base of SMFLG ||9|, [l^, 11, 12]. For a rather long fiber 



even small random inhomogeneities strongly affect radiation characteristics and make their 
theoretical analysis quite complicated. 

To investigate radiation characteristics in SMFLG with random inhomogeneities physi- 
cists usually apply methods of perturbation theory |1^, 14, 15, 16, 17|. These methods are, 



however, essentially restricted both in the magnitude of fluctuations and the length of the 
fiber. This forces one to interpretate the results of such calculations for sufficiently long 
fibers. As an illustration we point at the discussion about the value (zero or not zero?) 
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of the limit polarization degree in infinite fibers (see, e.g., [|T3|, U, |T5|, |l^). It was 

this very problem that initially prompted us to develop the methods used in this paper. 
In particular, we have proved that the limit polarization degree is zero (see Corollary |2.4| 
in Section 2). 

Here we propose a new approach to the mathematical description and investigation of 
light propagation in SMFLG with random inhomogeneities. Our methods provide with 
a correct analytical investigation of statistical characteristics of complex amplitudes of 
an electric field of light wave for arbitrarily long fibers. Note, nevertheless, that the 
mathematical model (equation (|2.1| )) we use is only applicable for real physical devices 
with not too strong fluctuations, because it does not take into account the reflection of 
light on inhomogeneities. The results, obtained in the framework of our model, present, 
however, the exact solution of mathematical problems formulated. 

Our model and methods give a convenient basis for numerical simulation of different 
optical fiber devises (see, e.g., [|^, |2^, 21, 22, p3|). 

Three pillars our approach is based on are: the random process theory, as an instru- 
ment for correct mathematical description of the statistics of inhomogeneities in fiber, the 
differential equation with random coefficients as mathematical model of wave propagation 
in random media and random groups theory as (and it seems to be quite unusual) an 
instrument for exact solving these equations. Of similar approaches to similar problems 
we only know [p^. 

On the structure of the paper. In Section 2 we introduce notions we need, formulate 
the problems we intend to solve and our main results. We also formulate two theorems 
concerning the random groups theory itself. Both these theorems are needed to solve our 
problems; one of them is only a reformulation in the form convenient to us of one of the 
main results of the random groups theory |2^], while the other one seems to be new, 
although very close to some known results from the Markov process theory |27]. Among 
"physically meaningful" results we state that the limit polarization degree is exactly zero. 
All the proofs, both of physical and pure mathematical nature, are collected in Section 3 
divided into subsections. Our calculations are sometimes quite complicated and sometimes 
simple but rather long so we separated the assertions from their proofs. We included 
more or less detailed proofs because we believe that our methods may be also of interest 
themselves. In proofs we omit some technical details when they do not contain in our 
opinion any nontrivial idea. 



2 Basic definitions, equations and main results 

We consider the propagation problem of a nonmonochromatic light wave in a single-mode 
fiber light-guide (SMFLG) with varying directions of anisotropic axes and without losses. 
Let z be the lengthwise coordinate in fiber, /? a parameter depending on wavelength 

and physical parameters of fiber, E{j3,z) = ( ^^j^^ | j the complex amplitude vector 

\Ey[fj, z)J 

of the two orthogonal components of the electric field of a monochromatic light wave at 
point z, where x- and y-axes coincide with directions of anisotropic axes in this point. 
The propagation of a monochromatic light wave in SMFLG is described by the following 
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equation for vector E{(3, z) ||T 
dE 



Xf3{e{z))E, Xf,{e{z)) 



V -e(.) -I 



(2.1) 



where 0(z) is the rotation velocity (twist) of anisotropic axes, i.e, the derivatives of the 
directions of anisotropic axes with respect to z. 

The real optical fibers have random inhomogeneities occasioned by uncontrollable pe- 
culiarities in the preparation process of fiber and fiber devices. We will assume in what 
follows that the function 0(z) is a realization of a random process. Observe that though in 
real fibers there appears not only fluctuation of anisotropic axes's directions, we disregard 
all other kinds of random inhomogeneities. The reason is that their influence on light 



propagation in fiber is not so essential as that of the considered one, see, e.g., 29, 3C]. 

Our main goal is to investigate the statistical characteristics of the solutions of equa- 
tion (^.ll) . To do so, we need a rigorous mathematics description of the process Q{z). It 
seems that the statistical characteristics we are interested in do not depend crucially on 
the details of this description. Therefore, we will consider a sufficiently simple mathema- 
tical model of such process, slightly more general than suggested in This model, in 
our mind, is close enough to description of the real structure of random inhomogeneities 
in fibers, and, on the other hand, allows us to obtain rigorous results. We give some addi- 
tional arguments in favour of our model in what follows. Note also that this model appears 



to be very useful in simulation of light propagation in random SMLFG [19, 20, 21, 22, 23]. 

Let {OfejfcLo {^A:}fcLo sequences of independent random values such that all 

Gfc are distributed with the same probability measure \xq on M and 1^ are also equivalently 
distributed with density pi supported on M^.. Let the random function 0(z) be defined 
by the relation 

fc-l k 

e(^) = efc, where ^lj<,z<^lj, and /c = l,2,.... (2.2) 

i=i i=i 

The random process ( p.2| ) is a simple and convenient mathematical model, its particular 
cases are often used in different applications of the probability theory [^]. As it was 
shown in [^] this model with specific [iq and pi is physically justified and results obtained 
in the frameworks of this model are supported by empirical observations. 

In what follows we denote by angular brackets ( • ) the mean value of a random value. 

We will say that the regular twist of the fiber is absent if the distribution of Q is 
symmetric with respect to zero, in particular, (0^"+^) = for all non-negative integers n. 

Let us now introduce parameters of the light field in fiber, the parameters which are 
most interesting from a physical point of view. 

The coherence matrix of a monochromatic light wave is defined by the relation 

J{P,z) = E{P,z)E^{P,z), (2.3) 



where f denotes the Hermitean conjugation. The mean coherence matrix {J{P, z)) is 
an important characteristic of the monochromatic wave, because the coupling (energy 
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exchange) of polarization modes that occurs in SMFLG with random inhomogeneities is 
described by the diagonal elements of this matrix. 

In what follows we will need the function B[[3) of spectral density of nonmonochromatic 
radiation in fiber by the relation 



B{f5) =tvJ{l3,z) 



(2.4) 



Note that due to equation ( ^.ID B{P) is indeed independent on z. 

A incoherent nonmonochromatic radiation in SMFLG is characterized by its polariza- 
tion degree p{z) p2|, which is 



p{z) 



4det J(z) 



tr2j(z) 



where J{z) is the coherence matrix of a nonmonochromatic radiation, 



J{z) 



J J{(3,z)d(3. 



(2.5) 



The mean square polarization degree is then equal to {p^[z)). The study of this value is 
one of the main goals of our work. 

Earlier, in the frameworks of perturbation theory it was shown (see, e.g., [^3|) that, in 

the absence of the regular twist and if ^(/3, 0) = ( ^ ) , the averaged intensities of the 



eigenmodes in SMFLG as z ^ oo are of the form 







where 



h{(3) = lim - 



{J22iP,z)) 



Q{z)e'^' dz 



-2h{/3)z 



). 



(2.6) 



(2.7) 



is the h-parameter used usually to characterize the coupling between polarization modes 
1^, 34|. As we will see in what follows, the relation like ( p.6| ) is valid in our model, 

but the formula for /i-parameter is slightly different from ( |2.7D and coincides with it only 
when the random twist is small enough. 

We wish to analyze the behavior of these characteristics both as z — > cxd and at fini- 
te z. It is easy to see that to calculate {J{f3,z)) and {p^{z)) it is sufficient to know the 
distribution of random vector E{P,z) and also joint distribution of the vectors E{Pi,z) 
and E{P2,z) at different Pi and /?2- We don't need other statistical parameters of vectors 
E{(3, z) and so we will investigate only the mentioned ones. 

Denote by N{z) the nonnegative integer such that 



N 



N{z) 



mm 



N 



> z 



(2i 
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Note that N{z), as a random integer function of z, depends on the reahzation of the 
random process Q{z). Now, the solution of equation {2A) with function Q[z) of the form 
( |2.2| ) can be expressed as 

E{P,z) = Uif3,z)Eoif3) 

I N{z)-1 \ 

= mAz- ^ ^fc,ejv(.)l M^(/^(,)_i,G^(,)_i)---M^(/i,Gi)^o(/3), (2.9) 

where Eq{P) is a complex amplitude of light field at the entry point of fiber, and {I, Q) 
= exp (/X^(G)) is Jones matrix of the fiber section of length / and fixed axis twist G, 



M^G,G) 



2G . flpg\ [ll3e\ .P . 

-—— sm cos — 1— sm , , 

V Pe \2 J \2 J Pe \2 J / 



(2.10) 



Here we set f3e = VP'^ + 4©^- 

Random matrices Mp{l, G) generate a subgroup G of the group SU{2) of 2 x 2 unitary 
matrices. Thus, the analysis of the statistics of E{I3, z) reduces to the analysis of statistics 
of the products of random matrices from some group, i.e., to a problem from the theory 
of random groups. 

In what follows we will need functions m^fc(/,G), k G {0, 1,3} such that 

M,(,e)= + ®' V (2.11) 

Together with the random vector E{f3, z) which describes the complex amplitude of the 
electric field at point z, we consider also a random vector i?jv(/3) which describes the field 
at the output of the fiber that consists of precisely random sections. The vector Ef^{f5) 
is defined by a formula similar to ( p.9D but where the matrix-valued function U{(5, z) is 
replaced with the matrix-valued function 

Un{I3) = Mp {In, On) Mp {l^-i, On^i) ■ ■ ■ Mp [h, Gi) . (2.12) 

We will see that it is possible to use this "discrete semigroup" instead of "continuous" one 
from (|2.9| ) in order to investigate the statistical properties of electric field in the random 
fiber. This is useful because there exist powerful results about the random matrix products 
in random group theory. 

Let us formulate now some useful theorems based on the random group theory. Let 
be a sequence of independent equivalently distributed on R+ random values and 
{cj}^i be a sequence of independent equivalently distributed random vectors with values 
in and absolutely continuous with respect to Lebesgue measure probability distribu- 
tions. Let also Y : M" — > g be smooth enough and almost everywhere (with respect to 
the measure fic generated by random vector c) nonconstant function on with values in 
the Lie algebra q of the orthogonal or unitary group acting in the m-dimensional space. 
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i.e., with values in the space of m x m skew-symmetric or skew-Hermitean matrices, re- 
spectively. Define the matrix-valued function A : M+xM" — > G, where G = SO{m) or 
G = SU (m), by setting 

A{z,q) = exp{zY{q)) . 

The matrix- valued function A{z,q) generates a random semigroup, i.e., the random 
matrix-valued function U{zi, Z2), where zi,Z2 G of the form 

/ N(Z2)-1 
U{z2,Zi) = aIz2- Y1 h^^N(z2) 




where N{z) is defined by ( |2.S| ) and Z2 ^ zi ^ 0. (Hereafter we assume Iq = 0.) The term 
"continuous random semigroup" is related with the fact that a solution of the stochastic 
differential equation ( p.lD is of such form, and also with the obvious fact that the equality 
U{z3, Z2)U{z2, zi) = U (^3, zi) holds for z^ ^ Z2 ^ zi ^ 0. Below we deal with the matrix- 
valued function U{z) = U{z,0). We note that the defined above random semigroup is 
stationary in the sense that the statistical characteristics of the matrix U{zi,Z2) depend 
only on Z2 — zi. 

Together with the continuous matrix semigroup, we will consider the discrete random 
semigroup Unk, N, K £ Z+ which is defined for integer N ^ K hy the equality 

UnK = A{In , Cn)A{In-1, CN~l) ■ ■ ■ A{lK+l,CK+l)A{lK,CK)■ 
As in the continuous case the term "semigroup" stems from the relation UnkUrl = Unl 
valid for N ^ K ^ L ^ the stationarity condition is realized and we will again be 
interested in the statistical properties of the family Un = Unq. 

We denote by G the minimal closed subgroup in G which contains all matrices {A{t, q) \ 
{t,q) £ M+x supple}- By definition, G is a non-discrete subgroup in the compact simple 
group G. 

This theorem is a specialization of "the central limit theorem" for compact stochastic 
groups for the considered case (see, for example, [|^, |2^). 

Theorem 2.1. Assume that the family of matrices {A{z,q) \ {z,q) £ R+x supp/^c} is not 
contained in any conjugacy class with respect to any normal subgroup of G. Then the 
probability distribution of the random matrix Un tends to the Haar's measure on G as 
N ooQ. 

The following theorem includes several statements which make it possible to calculate 
limit mean values of different linear (with respect to the semigroups introduced) quantities. 



Formally, this theorem does not depend on Theorem 2.1 but in fact its hypothesis are 



equivalent to the ones of Theorem 2.1 and its conclusions are corollaries of Theorem 2.1 



^ The condition concerning the conjugacy classes is present in the central limit theorem ^ for 
arbitrary compact stochastic group. This condition is certainly realized if G is simple: no normal subgroups. 
As we will see below, if G is semisimple this condition is equivalent to the statement that the family 
{A{z,q) I {z,q) £ R+ X supp /ic} does not have any common eigenvector corresponding to the eigenvalue 
different from 1. 



Random Groups in the Optical Waveguides Theory 



497 



Theorem 2.2. Let Pq be an orthogonal projection on a common eigenspace correspond- 
ing to the common zero eigenvalue of the family {Y(q) \ q £ supped- Suppose that the 
characteristic function pi{X), A G C, of measure pi, defined for Re A ^0 by the relation 

/•oo 

pi{X)= / pi{z)e-^'dz = {e-^'), (2.13) 
Jo 

is holomorphic in the closed left half-plane and has a meromorphic extension to the whole 
complex plane. Then 

1) the limit mean of matrices Um as N ^ oo exists and lim {Un) = -Po; 

2) the limit mean of matrices U{z) as z ^ oo exists and lim {U{z)) = Pq. 

z— »oo 

This theorem allows us to find some limit values using discrete semigroup Un instead 
of continuous U{z). As we see below, this is sometimes more simple to do. 

The following statements are obtained by applying the limit theorems on random matrix 
semigroups in order to derive statistical properties of a radiation in SMFLG when the fiber 
length tends to infinity. 

Corollary 2.1. The limit distribution of En{P) as N ^ oo is the uniform distribution 
on the three-dimensional sphere. 

Corollary 2.2. Let Pi / ±02- Then the vectors En{(3i) and E]y{P2) o,re independent as 
N ^ oo. 



Corollary 2.3. lim(J(/3,z)) 




Corollary 2.4. lim {p (z)) = 0. 

Z— >CXD 

These results deal with the limit properties of light waves in fibers. The result of 



Corollary 2^ is well known for fibers without regular twist (see, e.g., ||3^ ). Corollary |2.4| 
seems to be most meaningful for physicists. Note in this connection that the problem of 
determination of the limit polarization degree is widely discussed in the literature up to 
now, see, e.g., ||, |^, |l8|, |l| . 

It is very interesting also to find out the dependence of the mean coherence matrix 
and polarization degree on the length-wise coordinate z. Some particular results in this 
direction are listed bellow. We emphasize that all these results have an asymptotic nature, 
i.e., are valid approximately for sufficiently large z. 

First, some preliminary notices. 

As mentioned above, the problems connected with continuous semigroup U (z) are often 
essentially more complicated than the ones connected with Un- This is why in what follows 
we deal with Un instead of U{z) and intend now to explain some reasons for possibility 
to perform such replacement. In fact, we wish that there were no essential differences 
between the asymptotic behavior of {UnO (^(-^(O)C) as ^ oo. We did not try to 
prove or even to formulate this assertion in rigorous mathematical manner, but we base 
our wishes on various numerical experiments which confirmed this "fact" ||2^. We think 
that it is possible to make a theorem from this "intuitive" assertion, but its proof must 
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be quite difficult. The reader can anticipate this from presented in Section 3 proofs of 
asymptotic assertions for the discrete semigroup; they appear to be long enough even in 
this relatively simple case. 

Proposition 2.1. Suppose that in a fiber the regular twist is absent and 2{m'^^{l, ©)) < 1, 
where mpj is defined in ( \2.1]\) . Then the diagonal elements of the "mean coherence matrix" 
^ { f^\\ r,+ i^^r,^ AT f„,if;n nr-oi r,+n r^nn A^) A) o„^+h — N (l) and h'paramctcr of 



Un{P)J{(3,0)UIj{P)) at large N fulfill relation (g^ with z 
the form 



1 



8G2 



/32 + 46- 



■sin2 (^1^/32 + 482^^^ . (2.14) 



It is interesting to compare relation (2.14) with (2.7). To do so, we must obtain an 
explicit formula for the "classical" /i-parameter. The following assertion holds. 



Proposition 2.2. 




e{z)e"^' dz 



4(62) 



sm 



IP 



(2.15) 



To compare the different expressions for /i-parameter, i.e., ( |2.14| ) and ( ^.15 ), we con- 
sider the case when pi{z) = exp (^~^^- Suppose that supple C [—©max, ©max] with 
©max ^ y^/32 + (/)~2. Then it is easy to see that the ratio of expressions ( |2.14| ) and ( p. 15 ) 



46^ 



. Thus, the conditions of applicability of expression (E7 



is of the form 1 + ^ ^f^-2_|_^i 
appear to be quite rigidly restricted. 

Proposition 2.3. The asymptotic of the mean square of the polarization degree is of the 
form 



{P 



N 



where B{(3) 



tr Jjv 



_ 1 

~ 3 

g(/3) 
/ B(/3) d/3 




d(5 + 



1 



Ar3/2 



(2.16) 



is a normalized spectral function, B{[j) is defined in 



(2-4 ), /(/?) is a rational function of the averaged polynomials inm^j and their derivatives 
with respect to (3 (for details see proof in Section 3). In particular, if the regular twist is 
absent, then 



fee{o,i,3} 
dmpi(l,Q) 



' dmpj{l,@ y 
. dp 



+ 



32 



dp 



mi3o{l,Q) - mpi{l,Q) 



dmpQ{l,Q) 
dp 



Observe, that beside the power series on N, the full asymptotic expansion for (pj^) 
contains also a finite sum of terms of the form ae~"^. (One can see it from the proof 
of this Proposition.) Unfortunately, we cannot explicitly estimate the values a and a in 
these terms using our approach, hence, we cannot estimate the accuracy of asymptotic 
(|2l^ ) as function of N. 
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3 Proofs of basic assertions 



3.1 Proof of Theorem 2.2 



To prove heading 1), introduce matrix S = {A), the mean of the matrix family {A{z,q) \ 
{z,q) G M+x supple}, i.e., 

S = j Pi{z) j exp(zy(g)) d^ic{q)dz. 

Since the families {/j}^^ and {cj^JL^i are independent, we immediately see that 

{Un) = 5^. 

Obviously, SPq = Pq, therefore, Pq = Pq and lim S^Pq = Pq. 

Let us show now that lim {E — Pq) = 0, wherefrom the statement required. To 

N—>oo 

this end, consider an arbitrary nonzero vector x G C™ such that Pqx = 0. The unitarity of 
A implies that PqSx = 0, too. Indeed, if Ax = x, then x = A'^x and, therefore, Pq = S^Pq. 
Conjugating both sides of this equality we see that Pq = PqS. It is clear that absolute 
values of all eigenvalues of the matrix S do not exceed 1 since 

\\S\\ = \\{A(z,q))\\^{\\Aiz,q)\\)=l. 

Now we will need a lemma. 

Lemma 3.1. (i) If the absolute value of the eigenvalue r] of S is equal to 1, then r] = 1; 

(a) the vector e is an eigenvector with eigenvalue 1 of S if and only if it is a common 
eigenvector of the matrix family {A{z,q) \ {z,q) G M+x supped with eigenvalue 1. 

Proof. Let r] ^ 1 he the eigenvalue of the matrix S and Er^ be the corresponding norma- 
lized eigenvector. Then there exists go G supp Hc such that 

^exp(zy(go))e^/0, (3.1) 

because otherwise for any q £ supp /j^c we would have had 
d 







exp{zY{q))er, 



llexp izYiq))Yiq)eJ = \\Yiq)eJ, 
i.e., Y{q)£fi = and = 1 in contradiction with assumption. By continuity of Y the 



inequality (3.1) holds in some neighborhood ilq^ of positive /ic-nieasure of the point go- 
There are two possibilities: either is a common eigenvector of the matrix family 
{^(9) I Q £ supple} or there exists a go such that the vector exp (z y(go)) is non- 
collinear to £n for all nonzero z. 

Consider the first possibility. Since Y is skew-Hermitian or skew-orthogonal, all its 
eigenvalues are imaginary. Let ix{q), where xil) ^ be the eigenvalue of Y{q) corre- 
sponding to the eigenvector e^. From the condition xiQo) foi^ some go 'ind a continuity 
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of Y it follows that we can find a neighborhood Qqg so that the inequality |x('7)| > ^Ixili 
holds for any q E 0^^. Let pi{u!) be a characteristic function of the distribution pi, i.e., 



Pi 



It is known that |/0;(lj)| ^ 1 for all lo and < 1 for u; 7^ if is piecewise continuous 

function. So, \piixiQ))\ < 1 — i5 for all (7 G ^Ig^ and some positive 5. Then 



\r]\ = \{Ser„er,)\ = 
^ I \plixiq))\ dpciq) 



piixiq)) dpc{q) 



+ 



90 



\piixiQ))\ dp-c{q) ^ 1 - 5/ic( ^qo) < 1 



10/ 



as was required. 

The second possibility is considered in a similar way. Let zq > be such that Pi{zq) ^ 0. 
Then, by assumption, the vector exp (zq ^^(q'o)) is noncollinear to and for some 
5 > we have |(exp (zq Y^qq)) e^, e^)| = 1 — 25. It follows from the continuity of Y that 
|(exp {zY{q))er,, = 1 - (5 for ah {z,q) G [zq - A, zq + A] x Qg^. Thus, 



M = \{Se,j,en)\= J dpdq) J piit) {exp {zY{q)) Sr,, Er,) dz 

^ 1 - 6 fici^qo) f^i {[zo - A, zo + A]) < 1- 
This completes the proof of Lemma |3.1|. 



These estimates and the theorem on the Jordan form of matrices imply that there exists 
an < 1 such that the inequality 

\\Sx\\ ^ 

holds for all x such that Pqx = 0. Thus, inequality HS^xH ^ and S'-invariance of 

the projection Pq implies now that 

lim = Po. 



This completes the proof of the first heading of Theorem 2.2. 

Let us prove heading 2). The independence of the random values {cj}'jLi 
implies that 



{u{t)) EE v{t) =Y. w df^ciq,) L n Pl(^^ 

k=l-^ j=l •^E'j<Zj=i 



dL 



k-l 



k-l 



where 



l-OO 

$/(z) = / pi{s)ds. 

J z 



(3.2) 
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Introduce the matrix 

L{z) = I e^^('')d/.e(g), 
and calculate the Laplace transform V{X) of the matrix- valued function V{t), 

l-OO 

V{\) = / V{z)e-^'dz. 
Jo 

To do this, we use the convolution theorem for the Laplace transforms. Thus, setting 

pec poo 

Lp{X)= / pi{z)L{z)e-^' dz, L$(A)= / '^>i{z)L{z) e-^' dz, 



Jo 

we obtain the expression for ^(A) of the form 

k=l 



oo k—1 —1 

V{X) = (£p(A)) ~ i$(A) = [e- L,{\)) L$(A). 



Note that by the definition ||Lp(A)|| < 1 for Re A > 0; hence, all the singularities of the 
matrix- valued function in last formula are only poles lying in left half-plane and only A = 
is a pole with the zero imaginary part. (Recall that as follows from our hypotheses about 
density pi, the matrices Lp[\) and -L$(A) are meromorphic in the whole plane and regular 
in the closed right half-plane.) 

We now calculate the limit of V{z) as z — > oo via the inverse Laplace transformation. 
For a > we have 

1 pa+\oo 

lim V{z) = lim / V{\) e^' dX = Res V{X). 

z^oo z~*oo 27ri Ja^ioo -^=0 

To calculate the residue in this formula, observe that by hypotheses the matrix L{z) has 
a block structure, L{z) = Pq (B L-^{z); hence, Lp{X) = pi{X)Po © L^(A), where pi(X) is 
defined by Thus, 

V{X) = (1 - piiX))-^ PoL^iX) © V^{X). 

Only the first term in this expression has a pole at A = 0, and the corresponding residue 
is 

A=o ' '^"^ dpi{X) 

It is evident that 

dpijX) 

A=0 



Res V(X) 

A=0 

dpiiX) 



A=0 



dX 

and since PqL{z) = Pq, we obtain 



dX 

zpi{z) dz = (/), 



poo 

Po£$(0) = Po / M^)dz = {l)Po. 
Jo 



Thus, Res ^(A) = Pq which completes the proof. 

A=0 
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3.2 Proof of Corollary ^ 

The absolute value of vector En{P) does not depend on N due to (|2.lD . Thus, a distribu- 
tion of the vector Ej\f{f3) has support on the three-dimensional sphere of radius \Eq{P)\. 
The limit distribution of the vector E]\f{P) = U]\f{(3)EQ{f3) as N ^ oo for an arbitrary 
initial vector Eq{P), is determined uniquely by the limit distribution of the random matrix 
Un{P). 

The random matrices M^(/,0) defined by expression ( |2.1[1| ) are elements of the group 
SU{2). Let us show that the subgroup Ai C SU{2) generated by the matrix family 



{Ma(/, 0) I (Z, 0) € M+x supp/ie} coincides with SU{2). Then, by Theorem 2A, the limit 
distribution of Un{(3) defined by (|2.12| ), is the Haar measure on SU{2). The matrices 
depend on two independent random parameters: / and G. Thus, dim(A^) ^ 2. Since 
dim{SU{2)) = 3, either dim{M) = 2 or dim(A^) = 3. In the last case M = SU{2). It is 
known [^], that there is no two-dimensional subgroups in SU{2). We are done. 

Thus, the limit distribution of the random matrix U^iP) as — > oo is the Haar 
measure on SU{2); hence, the limit distribution of the vector En{P) = U]\f{P)EQ{P) is 
uniform on the three-dimensional sphere with radius \Eq{ 



3.3 Proof of Corollary ^ 

We will show that the limit joint distribution of the matrices Un{Pi) and UNih) is the 
Haar measure on the group SU{2) x SU{2). This means that these matrices and therefore 
the vectors E]\f{Pi) and Ei\f{P2) for /3i ±/?2 became statistically independent as — > oo. 

Let G be a subgroup of SU (2) x SU (2) . The elements of G are pairs of matrices {gi , 52) > 
51,2 G SU{2). Denote by Pr^ : G — > SU{2), i = 1,2, the homomorphic projections, 
i.e., Prj ((51, (72)) = di- The following statement is known from the theory of semisimple 
compact groups (see, for example, p6|]). 

Lemma 3.2. Any automorphisms of the group SU{2) is of the form Aut (5) = hf{g)h^, 
where h is some fixed element from SU{2) and f{g) is equal to either g or the complex 
conjugate matrix 'g. 

As a corollary we obtain 

Lemma 3.3. // the subgroup G C SU{2) x SU{2) is such that both Pri and V12 are 

epimorphisms, then either G = SU{2) x SU{2) or G = {((7, Aut (g)) \ g G SU{2)}, where 
Aut is an automorphism of the group SU{2). 

Proof. Since Prj(G) = 5*^7(2), then e G Prj(G) where e is the unit element of SU{2). 
Denote: Ge = ^^i^{e). Then Ge is a normal subgroup in SU (2) as the kernel of the 
homomorphism Pri. 

Set now Gj = {5 G SU{2) \ {e,g) S Ge} C SU{2). It is clearly that Gj is a subgroup; 
let us show that it is a normal one. Indeed, ghg~^ G SU{2) for any h € Ge C SU{2) and 
given g G SU{2) because SU{2) is a group. Since Pr2 is a surjection, there exists g such 
that (g, g) £ G and since Ge is normal, we see that 

{9, 9){e,h){g, g)~^ = {e,ghg'^) G Ge- 

Hence, ghg^^ G Ge for any g G SU{2), so Ge is a normal subgroup. 
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But SU{2) is a simple group, so either Gj = SU{2) or Gj = e. 

Consider the first case. Let gi S SU{2). Then since Pr2 is onto, there exists g' G SU{2) 
such that {gi,g') G G. Let now g2 € SU{2); then {e,g2{g')^^) € Gg C G, since 
92(5')"^ e '5C/(2) = Ge'. Therefore 

(6,52(9')"^) • (51,9) = (51,52) G G. 

This means that G = 5^7(2) x SU{2). 

Now, consider the case when Ge = e. In this case for each g £ SU (2) there is only one 
g' G SU{2) such that {g,g') £ G. Indeed, let there exist g" £ SU{2) such that g' / g" and 
{g,g") G G. Then {g,g')^^ • {g,g") G G since G is a group. On the other hand 

(5,5')^' • (9,9") = {glMg')-') ■ (91,9") = {e,igT'9") S G^. 

Therefore, {g')-^g" G Ge' = e and g' = g" . 

Hence, there exists a map / : SU{2) — > SU{2) such that f{g) = g' if and only if 
{9-: 9') £ G. The map / is a homomorphism. Indeed, let {g,g') G G and {h,h') G G. 
Then (gh,g'h') G G because of G is a group and by the definition of / we see that 
f{gh)=g'h' = fig)f{h). 

Let us show that / is an injective homomorphism, i.e., if gi 7^ 52, then f{gi) 7^ f{92)- 
The kernel of / is a normal subgroup of SU{2). This means that either ker/ = e and in 
this case / is an injective homomorphism, or ker / = SU{2). In the last case the image of 
SU{2) is the single element e which is impossible because Pr2(G) = SU{2). Since Pr2 is 
onto, / is a surjective homomorphism. So / is both an injective and surjective, hence, / is 
an automorphism and G = {(g, f{g)) \ g G SU{2)^. Proof of Lemma is completed. ■ 



Let now G C SU{2) x SU(2) be the group generated by pairs of matrices (Mg^ (/,©), 
Mp^{l,Q)). Thus, by Theorem |2.1| , the limit distribution of the matrix pair {Un{(3i), 
Un{P2)) as — > 00 is the Haar measure on G. Let us show that in our case G = SU{2) x 
SU (2). Since random matrices Mp{l, 0) generated the whole SU (2), then G satisfies to the 
hypothesis of Lemma |3.3| . Assume that G is different from SU{2) x SU{2). This means that 
there is an automorphism ^ of the group SU{2) such that G = {{g,'i'{g)) \ g G SU{2)}. 
In this case the Haar measure on G is the image of the Haar measure on SU (2) with 
respect to the natural isomorphism g {g,"^{g)) of groups SU{2) and G. In accordance 



with Lemma 3^, there exists a matrix h G SU{2) such that g2 = /i/(gi)/i^ for any pair 
{91,92) e G. 

Assume that f{g) = g. For a given complex matrix R we have 

giRgl) = (giRhg\h)) = (^J^^^^^gRhgUg^ h\ 

where (•) denotes the averaging over Haar measure on SU{2). Any element g G SU{2) 
may be uniquely represented in the form 

' r, h \ 

(3.3) 
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where a and h are complex numbers such that |ap + = 1. Introduce the angles ijj 
and from the relations 

a = e'*^cos??, 6 = e^^sinT9, i? G [0, 7r/2] , y?, G [0, 27r] ; 

and denote the elements of the matrix Rh by Vmn- Then 

, / aarii + abri2 + abr2i + bhr22 -abru + a^ri2 - 6^r2i + a&r22 \ 

^fi/l^r' = __ _ _2 — - - - - . 

y -abru - b ri2 + a r2i + a6r22 bbrn - abri2 - abr2i + aar22 ) 
We need now to calculate 

^Q/ci-fc2^fc3 5fc4^ = ^(cOST?)*^l+*=2(sint?)fc3+fc4gi{fcl-fc2Vgi(fc3-fc4)^^_ 

The invariant normalized measure on SU{2) is dg = sin 2^} (I'd dip di/j, therefore, 

if A;i 7^ k2 or ^ k^, 

a-aH'^b'A = { k,\k,\ (3.4) 



{ki + ^3 + 1) 



if ki = k2 and kj, = k^. 



Using ( p.4| ) we see that (^gRhg"^^ = ^ tr {Rh)E, where E is the unit matrix. Therefore, for 
R = we have 

[gih-'gl) = hi (3.5) 



Now we find the mean value of the product giRg^ in a different way. In the space 
Mat(2;C) of complex 2x2-matrices, consider a linear operator V of the form 

V{R) = Mp^{l,Q)RMl^{l,Q). (3.6) 



Then {Uj^{pi)RU]^{p2)) = (^) (^) and {giRglj = ^lim {V)^ {R). It is easy to see 
that V preserves the inner product in the space of complex matrices given by the formula 
{Ri,R2) = tr{RiR\). Therefore, y is a unitary operator. Consider the matrix- valued 
function R{1) = M^^ (/, G)i?Mj^ (/, 9). Since M^{l,e) = exp(/X^(e)), function R{1) sat- 
isfies the differential equation 

^ = X^_^ (G)i?(/) + R{1)XI (9) = Y{e)R{l) 



with the initial condition i?(0) = R. This means that operator (3.6) is of the form 
y = exp(^y(9)). 

The matrix (V) satisfies the hypothesis of Theorem |2.2| , so the limit of the expression 
{V)^{R) as — > oo exists and is equal to Po{R), where Pq is the projection on the common 
kernel of the operators 1^(9), or equivalently, the common eigenspace of the operator 
family (3.6), corresponding to the eigenvalue 1. Let us show that the operator family (3.6) 
does not have the common eigenvalue 1 when /3i ^ ±/32. It is obvious that if a matrix R 
is an eigenvector of operator ( |3.6|) with eigenvalue 1 for any pair (/,9) G M+x supple j 
then the identity 

Mfs,{l,e)RMl{l,@) = R (3.7) 
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holds. Differentiating both parts of 
= exp (I Xfj{@)), we see that 



with respect to I and using the relation Mf^{l, B) 



Xp,{Q)R + RXUQ) = Q. 



The matrix Xfj{Q) defined by formula (2.1) can be expressed as 



where ai 



and (72 



1 
-1 



CJl + 0(T2, 



. Then 



i/32 



0. 



Differentiating ( p.SD by we obtain 

(72 R- Ra2 

If we substitute ( |3.9| ) into ( |3.8| ) we see that 

PicriR - p2R(Ti = 0. 



(3.8) 



(3.9) 



(3.10) 



It follows from relations ( |3.9| ) and (3.10) that the elements rij of the matrix R satisfy the 
following system of equations 

rii - r22 = 0, 
ri2 + r2i = 0, 
iPi - P2)ru = 0, 
I (/?i + /32)ri2 = 0. 

It is obvious that this system has a nontrivial solution only if /3i = ±^2 and, therefore, 1 
is not a common eigenvalue of the operator family (3.6) when /?i ^ ±/32. 



Thus, it follows from heading 1) of Theorem 2^ that 



aiRgl 



hm {Un{Pi)RUJ,{P2: 

A'— >CXD \ 







for /?! 7^ ib/?2 and any complex matrix R in contradiction with ( |3.5| ). 

The case when the automorphism defining G contains complex conjugation can be 
considered in similar way. 

So the joint distribution of matrices UooiPi) and U^{j32) is the Haar measure on 
SU(2) X SU{2), i.e., matrices f7oo(/?i) and UooiP2) as well as vectors EodPi), Eoo{P2) 
are independent for different /3i / ±/32- 



3.4 Proofs of Corollaries and ^ 

Observe that, due to equation (^j|), i/ie tote/ energy of radiation, being equal to / B{P) dP, 
where B{P) is defined by ( |2.4D , is preserved during propagation of the radiation in SMFLG, 
so we may assume that / B{P) d(3 = 1. 
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Taking into account that E]y((3) = Un{I3)Eq{I3) we rewrite expression (2^) for the 



coherence matrix at the output of the fiber that consists of N random sections as 

Jn = IUj^{P)E^if})E^{(3)ul,{l3)d(3. (3.11) 

Since the trace of the product of two matrices does not depend on the order of factors, 
and the matrix U]\f{j3) is unitary, we deduce from ( |3.11 ) that 



tr = j tr (4t(/3)C/^(/?)C/^(/3)4(/3)) d(3 = J B{(3)d(3 = 1. (3.12) 

Thus, the trace of the coherence matrix at any point of the fiber is equal to 1, irrespectively 
of the concrete structure of inhomogeneities in the fiber. 

To calculate the limit mean value of the coherence matrix, consider the bilinear form 



{JNS,r)= / {s,En{(3)^ (r,ENi(3))df3, 



and find the limit mean value for this form for arbitrary vectors s and r from the space 
of two-dimensional complex vectors . Since the distribution of the vector En (/3) on the 
three-dimensional sphere of radius B{l3) as ^ oo is uniform (Corollary |2.l| ), we have 



{{Joos, f )) = J B{I3) dp J^^is, C){f, C) d^= f{s, f ), 

where 5"^ is the unit three-dimensional sphere, and / is a function of two vectors- arguments. 
It is easy to show that / is invariant with respect to an arbitrary rotation specified by 
the unitary matrix V, i.e., f{Vs,Vf) = f{s,f). Additionally, this function is linear in 
the first argument and antilinear in the second one. We will show now that /(s, r) = 
a(s, r), where a = const. Fixing the first argument of / and using linearity, we see that 
f{s,f) = ((7(s),r), where g{-) is a linear function on a finite dimensional vector space. 
Hence, there exists a matrix G such that g{s) = Gs. Thus, f{s,r) = {Gs,r). Using 
the invariance of / we have {GVs,Vr) = {V'^ GVs,r) = {Gs,r), i.e., G = V'^GV for any 
unitary matrix V . Due to the Schur lemma all eigenvalues of the matrix G with this 
property coincide, i.e., G = aE, where E is the unit matrix. Thus, ((JooS, r)) = a(s, r), 
where a is a constant. To find this constant, we calculate the limit mean value of the trace 
of the coherence matrix. Using the formula 

tr Jat = (JArei,ei) + (JAre2,e2), 

where {ei, 6*2} is an orthonormal basis in C^, we obtain: (tr Jqo) = 2a. But, on the other 
hand, since ( |3.12| ) we deduce that tr Jqo = tr Jq = 1. Therefore, a = ^ and 

{{Jo.s,r)) = l{s,f)- (3.13) 



It follows from ( |3.13| ) that ( Joo) = \ E. From the second part of Theorem 2J2 it follows 
also that lim {J{z)) = E. This is exactly the assertion of Corollary 2.3 . 
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Let us calculate now the limit mean value of the square of polarization degree. The 
square of the polarization degree at the output of the fiber that consists of random 
sections is described by the formula 

2 4det Jjv 

P^ = l- (3.14) 

Therefore, taking into account ( |3.12| ), to find the limit mean square of the degree of 
polarization using formula ( |3.14| ) , we have to find the limit mean value of det ■ To this 
end, we introduce an orthonormal basis {ei, 6*2} in C^. Then 

det Jat = (JArei,ei)(JAre2,e2) - |(JAre2, ei)|^. (3.15) 

Clearly, to calculate (det Jn), one should know how to find {{JNS,f){Jj\fU,v )) for arbitrary 
vectors s, f, u and v. We have 

{JNS,r){JNU,v) = (s,^jv(/3i)) (r,^7v(/3i)) (u,En{P2)) (vj^m) df3idP2- 

Since the vectors E]\j{/3i) and En{[32) are independent as ^ 00 (Corollary |2.2[ ), and 
their limit distributions are uniform (Corollary ^T| ), we obtain 

{{JooS,r){J^u,v)) = jj B{/3i)B{(32)d/3id/32 JJ ^{sX){K^Ku,t){^) d^dr 

. ' ' 1 



Bi(3) dpj J^^is, 0(r, d^jj^^ 0(^, ? ) dC= -^{s, f){u,v). 

Applying this result to relation ( p.l5| ) and taking into account that vectors ei and 6*2 are 
orthonormal, we see that (det Jqo) = \- Thus, it follows from ( |3.14|) that (p^) = 0, and 
since {poo) ^ a/ {pIo), we also have (poo) = 0. 

From the second part of Theorem 12. 21 it follows that lim {p'^(z)) = 0. 



3.5 Proof of Proposition |2.1 



We wish to estimate the dependence of the diagonal components of Jn{P) on A'' in the 
case when the distribution of @ is symmetric with respect to zero, i.e., (0'^'^"^^) = for all 
integers k. 

Clearly, the matrix Jn{P) = E^{P)Ej^{(5) is Hermitian. According to our model, the 
vector Ek at the endpoint of the A:-th fiber section results from the vector -E^^i at the 
endpoint of {k — l)-st fiber section by multiplying by the unitary matrix M^{li^, 0^) (see 
(|2l0|)). Hence, 



JkiP) = Mp{h,Qk)Jk-imMl{h,Qk) = M{(5- lk,ek){Jk~m)- (3.16) 
Now, consider the introduced linear operator in the space 7i of Hermitian 2x2-matrices. 



It is easy to see that the operator family (3.16) preserves the usual inner product (A, B) 



tr(A5t) in the space of complex matrices and, therefore, it is unitary. Observe also that 
operators from this family belong to the tensor product of the standard two-dimensional 
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representation of SU{2) in and its Hermitean conjugate; hence, have a unique common 
eigenvector '-'^o = ^ q ^ ^ with common eigenvalue 1. One can see also from ( 0.16 ) 
that if the distribution of is symmetric with respect to zero, the operator (M(/3)) has 
another, orthogonal to previous one, eigenvector o"i = ( ^ ? ) • Then the subspace of 



. -1 . 

the diagonal matrices is invariant with respect to the operator (M(/?)). Hence, we obtain 
the relation 

diag((J7v (/?))) = ^ (diag(Jo(/3)), ao) + ^ r?f (diag ( Jo(/?)), ai) ai, (3.17) 

where r]i is the eigenvalue of the operator (M(/3)) corresponding to the eigenvector ai. 
Clearly, 

r/i = i ((M(/5))ai,ai) = ^ (tr M^(/, GViM^ e)ai) = 1 - 2(m23(/,e)> . 
It follows now from (3.17) that 



{Mfi. N)) = 'I^(l- ,f ) . (3.18) 



tr Jo 

If rji > 0, this relation can be re-written in the form similar to (| 

Jo f , Jll — J 22 -2hz 



( J22(/?, ,)) = ^ f 1 - 'lll-lJk ) , (3.19) 



2 V tr Jo 

where z = N{1) and h = — 2^1nr/i. 
3.6 Proof of Proposition |2.2 



As in the proof of Theorem |2.2| , to calculate the value of the /i-parameter, we use the 

k 

Laplace transformation. Let N{z) be defined in ( |2.8D , z^ = X] h be defined in 

i=i 

(3.2). We deduce from ( |2.7] ), due to independence of random values lk,Qk, that 



/i(/?)=mn-/ ^ / / el^^e'f^^'~'Utds+ / e^(,)+ie'^(*-^)(itds 

Z—tCO Z \ J J J ZN(zy ^N(z) ' 



k 
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Using now the convolution theorem for the Laplace transformation, we can express the 
last expression as 



4(92) 1 



/32 
4(9^) 



27ri z 



lim — r / A 

;-*oo 27ri 



l-pi{\) {l-pi{\))\ 



/32 t^l\l-pi{\)^ {l-pi{X)f 



where /5;(A) is defined in ( 2.13 ) and 



OO -1 

^i{z)e-^' dz = -{l-pi{X)), 

^ 

OO 



Jo ^ 



pj(z)sin2^e-^^ dz. 



We see that ^z(A), V'i(A) and have finite limits as A ^ 0. Hence, 



/ AM/^A) AcDKA)V;KA) 
x=o [l - piiX) ^ (1 - piiXW ^ 

Thus, we obtain relation (|2.15|). 



MO) 



dpi{X) 



dX 
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3.7 Proof of Proposition pT3 



In this section we consider an asymptotic expansion of the mean square of the polarization 
degree at large A^. We can only find the principal term of this asymptotic. 

Note that simultaneously with the calculation of the asymptotic expansion of the mean 
square of the polarization degree we obtain another proof of the relation (p^) = 0. 

Our problem reduces to the asymptotic estimate of the mean value of the determinant of 
the coherence matrix. To do this, it suffices to calculate the mean value the tensor product 
Jn ® Jn, because the determinant of 2x2-matrix is a linear function of the elements of its 
tensor square. Further, since Jjy is an integral over the parameter f3 of matrices J]\f{P), 
the mentioned tensor product is the integral over the pair {^1,^2) of the tensor product 
JN{(ii) ® Jn{P2)- It follows from ( 3.16|) that there exists a linear operator M{Pi, (32] I, 9) 
m.n'^ = n®n such that 



Jn{I3i) ® Jn{(32) = M(/3i,/32; In,Qn) (JnMPi) ® JN-iiP2)) • 



(3.20) 



Note that the operator family M(/3i,/32; /,9) is unitary because it preserves the inner 
product {Ai(g> Bi,A2(g> B2) = {Ai, A2){Bi, B2) in H"^. 
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Hereafter we denote by {aj, j G {0,1,2,3}} the usual Pauli basis in the space H of 
Hermitian 2 x 2-matrices: 

10\ /10\ /01\ f i 





' 






0-2 = ( 






0-3 = f 



Wc also denote by sj = the corresponding orthonormal basis. Remember that the 
inner product in H is defined as {A,B) = tiAB. For the set {xa} of vectors in a given 
vector space we denote by C ({xa}) its linear envelope. 

Let us now list some facts from the representation theory of the group SU{2). 

Fl) Let TT be a representation of the group SU (2) in the space TC defined by the relation 
7r{g)A = gAg^ , A € Tl. Clearly, this representation is reducible and its decomposition into 
irreducible components is of the form H = Hq © i/a, where Hq = C{aQ) and = 

{{ox, 02,0^}). 

F2) Let TT®^ be a representation of SU{2) x SU{2) in such that 

T:®^{gug2){A®B) = T,{gi)A(^T,{g2)B. 

Then from the decomposition above it follows that 

n'^ = Ho® Ho® Ho® Hs® Ho® H3 = H^q © H^^ © iJ|o e i^fg 

and this is the decomposition of vr®^ into irreducible components. 

F3) Denote by tt^ the representation of SU{2) in H^, obtained by restricting tt®^ 
onto the image of the diagonal embedding of SU{2) into SU{2) x SU{2). Then the 
components Hqq, Hq^ and H^q are still irreducible, but the component H^^ can be further 
decomposed. Indeed, consider the natural action of the group Z2 = {0, 1} on n'^ by twist: 
1{A® B) = B I® A. Then H^^ = H^^g © H^^^ is the direct sum of the spaces of symmetric 
and skew-symmetric tensors under this action. The space H^^^ can be further decomposed 
into the direct sum 




^33s = ^ 

These are all irreducible components of vr^. 

Thus, the complete decomposition of tt^ in HP' consists of two one-dimensional com- 

3 

ponents, spanned by vectors Eq = so <Si sq and £1 = Si <8) Si, three three-dimensional 

1=1 

components in the spaces Hq^, H^q and -fffaa and one five-dimensional component in 

We will also need some relations between vectors and operators in Ti.'^. 

Rl) Let T = eo + £1- Then for A,B eU, one has {A ®B,T) = tiAB. 

R2) A = i (eo - ei). Then for AeU, one has {A © A, A) = det^. 

R3) Let 5 : M — SU{2) be a smooth function, and 52 : — SU{2) x SU{2) be such 
that g2{x,y) = {g{x), g{y)). Let further D be a differential operator in C°°(M^) of the 
form D = ^ — (Here we denote by C°°{M?) the space of smooth functions with values 
in an arbitrary finite dimensional vector space.) The operator Do : C°°(M^) — > C°°(M) is 
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then D\y^^, and similarly = D'^\y^^. The following relations can be proven by direct, 
but long, calculations: 

Dl {^^\g2{x,y))e,, e,) = -4tr ^ ^, (3.21) 

^o^^'(52(x,y))eiGF|3,. (3.22) 

Denote by L{j3i,(32) the mean operator {M{Pi, P2', ^0)) and by -F/v(/3i, /92) the mean 
value of Jn{Pi) ® Jn{P2)- Then it follows from ( |3.20| ) that 

FN{f3i,/32) = L^{Pi,P2)Fo{Pi,l32), (3.23) 

and we can use, as in the previous section, the spectral decomposition of the operator 
L{Pi,P2) to calculate F]\f{Pi, (32)- 

Note that we deal here with operators in the 16-dimensional space TC^ and there are no 
reasonable reasons for existence of an analytical solution of spectral problem for operator 
L(/3i,/32). This is the main obstacle for obtaining complete asymptotic decomposition 
of the polarization degree. We will see, nevertheless, that it is possible to obtain some 
analytical expressions for two major eigenvalues and corresponding eigenvectors of this 
operator which suffices to construct the leading term of the asymptotic of the polarization 
degree. 

Using definitions above, we deduce: 

(det Jtv) = (F,v(/3i,/?2),A) dPid(32 

(l^(/?i,/32)Fo(/3i,/32),a) dMp2. (3.24) 



The integrand in ( 3.24 ) can be written in the form 



L^(/3i,/32)Fo(/3i,/32),a) = - [Fo{Pi, P2), (l^ {Pi, P2)) {eo - ei] 

= ^{Fo{pi,P2),eo)-^(^Fo{Pi,P2),{LHpi,P2)Y ei^ . (3.25) 

As follows from the definition of vector eq, 

(Fo(/3i,/32),eo) = ^tr Jo(/3i)tr Jo(/?2) = ^B{P,)B{P2). (3.26) 

Let {sliPi, P2), ViiPi, l32)}j=o be the set of normalized eigenvectors and corresponding 
eigenvalues of L^{Pi,P2). Then the dual basis is the set {ei(/3i, /?2)}|£o eigenvectors of 
L{Pi,P2). Moreover, it follows from previous discussions that the vector eo(/3i,/52) = 
Sq{Pi,P2) = eo is orthogonal to all other eigenvectors Sj, f?o(/3i,/92) = 1, ^\{PiP) = 
ei(/3, P) = and ?7i(/3, P) = 1. Further, it follows from Theorem ^ that |?/i(/3i, /32)| ^ 

V < 1 for i = 2, . . . , 15. Expanding Ei in terms of the eigenvectors e\ we obtain 

i=0 

15 



i,/32)) £1 = ^(el,e^(/?l,/32)) p(/3i,/32)) 

i=0 

15 

= E (/3i, /92) (ei, e^(/?i, /32)) /92). 



i=l 
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It follows from the properties of eigenvalues r/j that the relation 



N 



-^f^(P^'P^) (ei,£i(/3i,/32)) 4(/?i,/92) +0 (e-^) 



holds for a = — Ini/ > and /i(/3i,/32) = — lnr/i(/3i, /32). 
Substituting (|3.27|) into (p.25|) and using (|3.26|) we obtain 



1 



_ lg-Affe(/3i,/32) 



ei,ei 



Now by substituting (|3.28|) into ( p.24 ), changing variables 
and integrating over 6 by the saddle-point method, we obtain 

2 



(det Jn) 



1 



B{p)df3 



27r 



dp + 



Ar3/2 



where /(/?) = -Do^(/3, /?). Here we use the relation 

Doh{p,(3) = -DoVi{P,P) = -1^0 {L\p,p)ei, ei) = 0, 
which follows from the standard perturbation theory |37] and the fact that 



2j£l, £1 



£^(/32,/3i: 



ei, El 



(To prove the last relation it suffices to note that, by definition, 

(L^{Pi,f32){A<S)A), A^A^ 
[Up, {I, @)AmI (/, 9) Mf3, e)AMl^ {l,e),A(^A 



tr i^Mp,{l,e)AMl^{l,e)A) tr (^M^,(/, G)ylMj^(/, 9)^ 

for any Hermitian 2 x 2-matrix A.) Taking into account the relations / B{P) df3 = It Jn 
and El = £q — 2A, we can express ( |3.2g| ) in the form 



(3.27) 



(3.28) 



(3.29) 



(3.30) 



^ i?JN - (det JN) = lj Q tr2 Jo(/5) - det Jq 



2lT 



Nf{(3) 



d(3 + 



Ar3/2 



Using now the formula (3.14) and a normalized spectral function 

tr Jo(/?) _ B{p) 



trJiv jB{l3)d(3' 
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we obtain expression ( |2.16| ) for the polarization degree. 

Our last problem is to calculate the value f{(5) in terms of parameters of Q{z). To do 
this, observe that 

Dlr,^{P,(]) = Die-^^f^'''^ = -Dih{P,(3), 
because DQh{(5, (3) = 0. Let now apply to the relation 

i^(/3i,/52)ei(/3i,/32) = ??i(/?i,/32)ei(/3i,/52) (3-31) 
and then scalar it by £i. We obtain: 

/(/?) = Z)2M/3,/3) 

= -1 [DlL\(3,l3)ei, ei) - (z)o^^(/3, /3)Doei(/3, /?), ei) . (3.32) 
The first term in the right hand side of ( p. 321) is calculated by using relation (3.21), 

[DlL\(3,(3)ei, El) = {Dl (7r®2(M^(/, G), M^(/, G))ei, ei)> 

/ mJ(/,0) M^(/,G)\ , , 

To calculate the second term, we need to find the vector DQei{(3, (3). Applying Dq to 
equation ( 3.31 ) and taking into account relation ( |3.30| ), we obtain for this vector the 
equation 

(^E-L^{P,(]))Doei{p,l3) = ^DoLH(],(3)e,. (3.34) 

Due to relation ( p. 22 ), the right hand side of equation ( 3.34| ) belongs to the space -fffsa, 
which is invariant, as follows from F3), under the operator in the left hand side. The space 
i/lg^ is, by definition, orthogonal to the spaces H^Q and (see F2) and F3)). The direct 
sum of these spaces is exactly the kernel of the operator E — L^{/3, (3), hence, the equation 
(|3|) has a unique solution which is orthogonal to H^q © iJ^, and this solution belongs to 
i?|3„. Since vector ei(/9i,/32) is normalized and orthogonal to Eq, this solution is the one 
we need. So, we can rewrite equation (|3.34| ) as an equation in the three-dimensional space 
H^^^- It is easy to see that the orthonormalized basis in H^^^ can be chosen in the form 

Si (g) S2 - S2 Sl , S2 S3 - •SS "X" ■S2 , S3 (gi Si - Si S3 

= 71 ' = 7=2 ' = 7=2 • 

Introduce the 3 x 3-matrix Sij = (^(^E — {l3 , l3)^ (l)j,(pi^, i,j £ {1,2,3}, and vectors 
Vi = {DoEi{p,l3),4>i) and Ui = (DoL^f3, P)Ei,4>i) , i G {1,2,3}. Then equation ( pl 



takes the form 

Sv = ^u. (3.35) 
v3 



By direct calculation we see that 





"ll3 

"1-33 "Zo3 I , u 

-mo3 1 - 




514 



G Malykin, V Pozdnyakova and I Shereshevskii 



Here 



mkj = {mf3kil,@)mpj{l,@)), 
"■^j = \ ^fl rni3j[l, 9) - m/3fc(/, 9) 



d(3 



d(3 



m/3k{l,Q) = -trM^(/,9)afc, k,j G {0,1,3}. 

We can now express the function /(/?) in terms of components of vector v and elements 
of matrix M^(Z,9) and its derivatives with respect to (3. Using ( |3.32| ), ( |3.33 ) we finally 
obtain 



/(/?) 



E 

fce{o,i,3} 



dmfjj{l,e) 
d(3 



-vid2.Q + V2diQ - vzd2,i) 



This relation is essentially simplified when the regular twist is absent: u and many elements 
of S vanish and the solution of equation ( ^.35 ) attains a simple explicit form: 



\ 



-2 /- "^'^ 



3 




So, in this case 



/(/?) = ^ 



3 \ 4(/32 + 492; 
/ //32 



+ 



492sin2-V/32 + 492' 
(/32 + 492)2 



4 \/32 + 492 
^ 3 7 9^ 



492 sin/v^/32 + 492 Y 
+ (/32 + 492)3/2 I 

I 



This completes the proof of Proposition 
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